Reexamination of a multisetting Bell inequality for qudits 
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The class of d-setting, d-outcome Bell inequalities proposed by Ji and collaborators [Phys. Rev. 
A 78, 052103] are reexamined. For every positive integer d > 2, we show that the corresponding 
non-trivial Bell inequality for probabilities provides the maximum classical winning probability of 
the Clauser-Horne-Shimony-Holt-like game with d inputs and d outputs. We also demonstrate that 
the general classical upper bounds given by Ji et al. axe underestimated, which invalidates many 
of the corresponding correlation inequalities presented thereof. We remedy this problem, partially, 
by providing the actual classical upper bound for d < 13 (including non-prime values of d). We 
further determine that for prime value d in this range, most of these probability and correlation 
inequalities are tight, i.e., facet-inducing for the respective classical correlation polytope. Stronger 
lower and upper bounds on the quantum violation of these inequalities are obtained. In particular, 
we prove that once the probability inequalities are given, their correlation counterparts given by Ji 
and co-workers are no longer relevant in terms of detecting the entanglement of a quantum state. 

PACS numbers: 03.65.Ud, 03.67.Mn 



I. INTRODUCTION 

Bell inequalities ^, being constraints that have 
to be satisfied by classical correlations, have long played 
an important role in shaping our current world view 
With the advent of quantum information science, these 
inequalities have also found applications in the design 
of quantum key distribution protocol jj] and its security 
analysis Q, as well as the reduction of communication 
complexity More recently, there is also a growing 
interest in thinking about these inequalities in the form of 
non-local games ![ ^ which, in turn, are closely related 
to the studies of interactive proof systems in computer 
science (see, for example, Refs. Q and ||]). 

To date, the studies of Bell inequalities have focused 
predominantly on those involving only binary outcomes, 
such as the Bell-Clauser-Horne-Shimony-Holt inequal- 
ity |10, ll| and the Bell-Clauscr Hornc (henceforth abbre- 
viated, respectively, as Bell-CHSH and BcU-CH) inequal- 
ity ^ |I| (see, for example, Refs. (l5[ |l^ and 
references therein for a review on bipartite two-outcome 
inequalities). This is, of course, by no means accidental 
as many of the quantum information processing proto- 
cols have been developed explicitly with qubits, i.e., two- 
level quantum systems in mind |20[ . However, given that 
higher-dimensional quantum systems are gaining impor- 
tance in quantum information processing tasks [ pT| |2^ , 
the time is now ripe to also perform further studies on 
multiple-outcome Bell inequalities, which are naturally 
suited for higher-dimensional quantum systems. 

In this regard, we note that there arc only relatively 
few works devoted to the studies of such Bell inequalities 
and their quantum-mechanical violations. For an exper- 
imental scenario involving only two subsystems, the pio- 
neering work by Collins et al. |G3| resulted in a class of 



Bell inequalities that involves two multiple-outcome mea- 
surements per site (see also Refs. and [^). This class 
of inequalities, now known as the CoUins-Gisin-Linden- 
Massar-Popcscu (CGLMP) inequalities, is tight i.e., 
they represent boundaries of the corresponding set of 
classical correlations, or more precisely, facets of the re- 
spective correlation polytope ||27| (for a review on the 
subject of polytope, see Ref. [^8|). 

Apart from the CGLMP inequalities, there are only 
a few other classes of Bell inequalities that are specifi- 
cally catered for multiple outcomes. Some of these are 
defined in terms of joint and marginal probabilities of ex- 
perimental outcomes | [l 5 , | 2 ^ . |3C| , [3l| , 32, ^ whereas 
the others 36, |37|, 38|, |39|] arc defined in terms of 
correlation functions — i.e., expectation value of the 
product of experimental outcomes. In general, however, 
very little is known about the tightness of these inequal- 
ities @||||,||,|0|. 

An interesting feature of tight multiple-outcome 
Bell inequalities is that, except for those introduced in 
Ref. Jssf , they are typically violated maximally by non- 
maximally entangled states [ 
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, ||. This and 

other evidences gathered from the studies of non-local 
apparatuses [Q — collectively known as "an anomaly 
of non-locality" — have led to the proposal of see- 
ing quantum entanglement and Bell inequality violation 
as fundamentally different resources p6| , even though we 
now know that all bipartite entangled states cannot be 
simulated by classical correlations alone p7[ . 

In this paper, we reexamine the class of bipartite d- 
setting, d-outcome Bell correlation inequalities proposed 
by Ji et al. |8|. In Sec. ||, we rewrite these correlation 
inequalities as Bell inequalities for probabilities and show 
that it admits a natural interpretation within the frame- 
work of the so-called CHSH game (but now with d 
inputs and d outputs). In the same section, we provide, 
for d < 13, the actual classical upper bound and for the 
more complicated scenarios, some non-trivial estimates 
thereof. The tightness of these inequalities is discussed 
in Section III. After that, in Sec. IV, we investigate the 
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quantum violation of the probabilities inequalities and 
compare them against those obtained in Ref. using 
their correlation counterparts. We will conclude with 
a summary of results and some possibilities for future 
research in Sec. Throughout, our discussion focuses 
on the scenarios where d is a prime number; the analo- 
gous computational results for non-prime value of d (with 
d < 12) are summarized briefly in Appendix |^. 



II. THE BELL FUNCTIONS AND THEIR 
CLASSICAL BOUNDS 

The Bell inequalities proposed by Ji and co- 
workers [^8| are applicable to an experimental scenario 
where two spatially separated experimenters (hereafter 
called Alice and Bob) are each allowed to perform d alter- 
native measurements, with d being arbitrary prime num- 
ber. Specifically, if we denote by w = e '^^r/d ^^.^j^ 
root of unity, Ji et al. consider local observables As^ and 
Bg^ that are unitary so that each measurement admits 
the d possible outcomes In these notations, the 

Bell function and the correlation inequalities presented in 
Ref. Q — up to a factor of l/((i — 1) — read as 

'^■" = E E ^^'^'HiAsSiBsS), (la) 

71—1 Sa-,Sb—0 

d ALf„ < 5„ + d^ <d Ai±, (lb) 

where the classical upper (lower) bound is determined by 
maximizing (minimizing) over all extremal (determinis- 
tic) classical strategies Oa{sa) and Ob(sh), i.e., 

^ill = max 6s^,^+o^(^s^)+oUsi), (Ic) 

and ^l^l^mm^ 6s^s,+o^{s^)+o,{s,), (Id) 

and here, Sj is a shorthand for the Kronecker delta 

f^O, j mod d- 

In Eq. (|l^), ((As^)" (BsJ") is a correla- 
tion function that gives the statistical average of the 
product of measurement outcomes of {Ag^)" and (B^J". 
Hereafter, we shall refer to the inequality upper (lower) 
bounding Sj^ in Eq. (|^) as 

For any given setup of the Bell experiment, the set 
of correlation functions {((^s„)" (^si,)")}s_^ „ can be 
collected together and written as the entries of a vector 
in C'' It is known that the set of such vectors al- 

lowed by a local hidden-variable theory (LHVT) forms a 
convex polytope pTf — i.e., loosely, higher-dimensional 
generalizations of convex polygons — called a classical 
correlation polytope. Each I^^ given in Eq. (|l|) then de- 
fines a hyperplane in the space of complex correlations 
(^d separating (some) correlations not attainable 

using LHVT from the classical correlation polytope. 

Now, let us further denote by p^g''(sa,Sh) the joint 
probability of Alice observing the Oa-th outcome and Bob 
observing the Ob-th outcome conditioned on her measur- 
ing As^ and him measuring Bs^^ ; likewise for the marginal 
probabilities p'^{sa) and p°^{sb)- From Eq. (||) and the 



fact that classically, 

d-l d-l 

{{AsSiBsS) = E E t^"°°+"°«^(sa,s.), (2) 

o„=0 Oi,=0 

it can be shown that the following Bell function, 

^ d-l 

^ ^ 'cP E ^s^s,+o^+o,P°^B''{Sa,Sb), (3a) 

must also be bounded from below and above as follows: 

/,-:^aL11<5, /+:5<lAil. (3b) 

This gives rise to two classes of linear Bell inequalities for 
probabilities. 

A few remarks are now in order. In contrast with 
^td given in Eq. (|l|) , the Bell inequalities given in Eq. (^) 

live in the space of real correlations R'** where each vector 
in the space has entries given by all the distinct joint 
probabilities p'^^'' {sa, Sb). Moreover, it is also easy to 
see that the requirement that each experimenter must 
perform unitary measurements in Eq. (|l|) is now lifted 
from I^; Alice and Bob are free to assign arbitrary values 
to their measurement outcomes. 

On the other hand, note that S only involves non- 
negative combination of p'^^'' {sa, Sb) and that the right- 
hand-side of Eq. (^^ is upper bounded |Q by 1. As a 
result, S can also be seen as the win ning probability of 
a two-prover, one-round unique game [p9[ whereby Alice 
and Bob win if and only if the answers that they provide 
Oa{sa) and Ob{sb) for the questions Sa, Sb (supplied to 
them with uniform probability) are such that: 

SaSb + Oaisa) + Ob{sb) mod d = 0. (4) 

This is clearly a direct generalization of the CHSH game 
presented in Ref. . Classically, the winning probability 
of the CHSH game [corresponding to d = 2 in Eq. (^)] is 
upper bounded by 3/4, but one can easily check that this 
is just the requirement of the Bell-CH inequalities p2| . 
For the rest of the paper, we will thus focus on scenarios 
where d > 2 and analyze the probability inequalities 
in connection with their correlation counterpart I^^. 

What are the actual classical bounds for these Bell 
inequalities? Here, we follow Ref. ||3^ and consider a 
d X d matrix Aid with its (s^ + 1, Sf, + 1) matrix element 
given by the left-hand-side of Eq. (^. A given extremal 
classical strategy, i.e., one that satisfies, 

Pab" ^b) = Pa {sa)p°B (sb), (5a) 
p°j{sa)^ 0,1, pl'isb)^0,l, (5b) 

then gives rise to a classical value of S and 5,1 determined 
by the number of zero entries in the corresponding matrix 
Aid- For d > 2, the following classical strategy |Q 

OaiSa) ^ Sa - l\f Sa ^ 0, 0^(0) = 0; 

Ob{sb) = iy Sb<d-2, Ob(d-l) = 2 

gives 5 = 0. This, together with the non-negativity of S 
[cf. Eq. (H)] show that Aif„ = 0. Thus, and J-^ for 
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d > 2 are Bell inequalities that are trivially satisfied by 
any theories that respect the non-negativity of probabil- 
ities. 

As for the classical upper bound, it was estimated 
in Rcf. |§ to be A'^l = 3{d - 1). While their exphcit 
analysis for d = 3 is valid, it can be verified that the 
following classical strategy for prime value d > 5, 

OaiSa) = d - 1 V Sa 7^ 0, 1, ^ , 
Oa(0)=0, Oa{l)=d-i, Oa (^^-i-^^ = - 3, 

Ob{sb) = V S6 > 2, 
Ob(0) = l, Ofc(l) = 3, Ofc(2) = 2 (6) 

gives 3d — 2 zero entries in A^^, indicating that A[^}^ > 
3d - 2. 

In this regard, we note that the actual value of A^f^: 
for d < 13 can be determined by exhaustively searching 
through all (extremal) classical strategies with the help 
from the following observations (all arithmetic operations 
described below are to be evaluated modulo d): 

1. 5 is invariant under the simultaneous transforma- 
tions: OaiSa) OaiSa) + k, Ob{sb) — > Oh(sh) - k for 

all Sa, Sb and all k £ {0,1, . . . ,d— 1}. Without loss 
of generality, we may thus set Oa(0) = 0. 

2. For all k e {1,2, . . . ,d ~ 1}, the strategics 

{oa{sa),Ob{sb)y-X=o and {o'M^ o'^{ksb)y-X=o 
give the same S if o'^{sa) = koa{sa) and o'^^{ksb) = 
kob{sb). This follows from Sa{ksb) + ojj(sa) + 
o'^{ksb) = k[saSb + Oa{sa) + Ob{sb)] and thus the 
two strategies give the same number of zeros in the 
corresponding dx d matrix Add- As a result, it suf- 
fices to consider Oa(l) = and Oa(l) = 1 once we 
have set Oa{0) = 0. 

3. For a given Sb and a given choice of {oa{sa)}'^~lQ, 
let k be the number in {0, 1, . . . , d — 1} that occurs 
most frequently in the expression "saSb -I- Oa(sa) 
mod d" as Sa varies from to d — 1. The optimum 
choice of Ob{sb) for the given Sb is Ob{sb) = d — k. 

Explicit value of these A^^}^ can be found in Table |. 

TABLE I: Aifi, and its lower bounds. The first row gives 
the values of the parameter d. The second row gives the best 
lower bound on Af^^x that we have found whereas its actual 
value, if known, is included in the third row of the table. 
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5 


7 


11 


13 


17 


19 


23 


29 


31 




6 
6 


12 
12 


19 
19 


37 
37 


47 
47 


66 


79 


99 


135 


148 



a random choice of classical strategy for Alice, the opti- 
mal strategy for Bob can be decided using the third ob- 
servation mentioned above. With this optimized classical 
strategy for Bob, we can in turn determine the optimal 
classical strategy for Alice and so on and so forth. The 
explicit values for some of these lower bounds, which we 
will denote by A[5 can be found in Table |. 



III. TIGHTNESS OF BELL INEQUALITIES 

A natural question that follows is whether the Bell 
inequalities and Jj" are tight, or so called facet- 
inducing for the respective set of classical correla- 
tions. By analyzing the tightness of these inequalities, 
we can gain insights on the structure of the correspond- 
ing set of classical correlations (Fig. |^). To this end, 
we note that the relevant classical correlation polytope 
for resides in a subspace of M'' of dimension 
dp = d^(d — 1)^ -I- 2d(d — 1). A Bell inequality is facet- 
inducing if and only if the number of linearly indepen- 
dent extremal classical probability (correlation) vectors 
saturating the inequality equals to the dimension of the 
polytope. For , this can be shown to be dp following a 
similar argument as that presented in Rcf. pq |; likewise 
for which can be shown to be dc — d^{d — 1). 




FIG. 1: (Color online) Schematic diagram of a two- 
dimensional plane in the space of (quantum) correlations. The 
shaded (light blue) polygon only consists of classical correla- 
tions whereas the convex region marked by a circumscribing 
solid curve also consists of nonclassical correlations. Ii cor- 
responds to a trivial Bell inequality that cannot be violated 
by quantum mechanics. The analog of a Ught Bell inequal- 
ity, such as I2 is the hyperplane given by I2 (red solid line) 
whereas the analog of a non-tight Bell inequality such as 
is given by I3 (blue dashed line). A quantum correlation "x" 
violates a Bell inequality if and only if the corresponding hy- 
perplane (eg. I2) separates "x" from the set of classical cor- 
relations. 



For d > 13, it seems formidable to search through all 
inequivalent (extremal) classical strategies 50 ; neither is 
the classical strategy given in Eq. ^ optimal. However, 
non-trivial lower bounds on A^^}^ can be obtained by op- 
timizing the classical strategies of Alice {oa(sa)}f =0 and 
Bob {ob{sb)}'^~2o iteratively. Specifically, if we start with 



In this regard, we note that our investigation shows 
that (Table ||) for d < 13, most of these probability and 
correlation inequalities are indeed facet-inducing. Note, 
however, that there is a priori no reason to expect that 
the Bell function given by Eq. (|l|) or Eq. (||) would give 
rise to any tight Bell inequalities. 
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TABLE II: Computational results for the tightness of and 
Jj" for d = 3, 5, 7, 11 and 13. The first column of the table gives 
the parameter d. From the second to the fourth columns, we 
have, respectively, the dimension of the correlation polytope 
relevant to J^t^j, the number of linearly independent extremal 
classical correlation vectors saturating inequahty J^^, and the 
tightness of The analogous results for jj", when available, 
are presented from the fifth to the seventh columns. 
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18 
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non-tight 


48 


18 


non-tight 
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100 


100 


tight 


440 


440 


tight 
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294 


294 


tight 


1848 


1848 


tight 


11 


1210 


1210 


tight 


12320 


12320 


tight 


13 


2028 


2028 


tight 


24648 







IV. QUANTUM VIOLATIONS 

In this section, wc will investigate the quantum vio- 
lations of /j" and compare them against those presented 
in Ref. |3|]. These quantities put bounds on the set of 
quantum correlations . In particular, the maximal vi- 
olation of a tight Bell inequality for a given state p is a 
primitive measure of the extent to which p is nonclassical. 
For example, in Fig. [l], if a Bell inequality is such that 
its Bell function gives zero for all points lying on I2 , then 
the maximal extent to which p violates this inequality, de- 
noted by 5qS-*(p), indicates the largest possible distance 
between an X2-violating correlation derivable from p and 
the hyperplane X2. Likewise, the largest possible distance 
between any point on the arc opposite to the polygon and 
I2 gives rise to the maximal possible quantum violation 
of the Bell inequality corresponding to T2 ■ 

Now, let us start by comparing the strength of 
against in terms of detecting nonclassical correlations 
present in an entangled state. To this end, it is worth not- 
ing that if we denote by the positive-operator-valued 
measure (POVM) element associated with the Oa-th out- 
come of Alice's Sa-th measurement, the expression 



d-l 

Oa.=0 



d-l 

E 



(7) 



holds true if and only if all the POVM elements satisfy 

A°'^A^Sa = ^Oao' A""^. This implies, in particular, that in 
quantum mechanics, Eq. (|^) is only applicable when we 
are considering projective measurements. In this case, it 
is easy to show using Born's rule, Eqs. ( |la| ) and (|3a|), that 
for any quantum state p, the quantum values of their Bell 
functions are related by 



tr(pS„) = d^tr(p6)-d2 



(8) 



where By, and B are, respectively, the Bell operator ||53| 
constructed from the Bell inequalities and . For 
generalized measurements where Eq. does not hold, 
measuring As^ no longer measures (A^^)" concurrently; 
Bji is also generally non-Hermitian in this case. Clearly, 
this makes a test of the quantum mechanical prediction 



against its classical counterpart [cf. Eq. (|l^)] meaning- 
less Thus, any quantum state that violates the in- 
equalities presented in Ref. |Q (and hence I'^j) must also 
violate /j" but the converse is not necessarily true. 

Numerically, by maximizing over the set of rank- 
1 projective measurements realizable through symmet- 
ric multipart beam splitters (see Ref. 55 and refer- 
ences therein), we have obtained some lower bounds 
on the maximal violation of with c? < 13 for 
the c?-dimensional maximally entangled state l^'J) = 
TsEfcTo where |u) and lig) are, respectively, 

the «-th computational basis state of AHce and Bob's sub- 
system. Explicit value of these quantum-mechanical vi- 
olation for \^'^), which we will denote by Sq^]^^{\'^'^)), 
can be found in Table III. Note that except for d = 5, 



these values also represent the best violation by \^~^) 
that we were able to find. Moreover, for d = 3 and 7, 

iSqM.'mpd^'jl)) are in fact the largest quantum violation 
of It that we have found. 



TABLE III: Bounds on the maximal quantum violation of /j". 
The first two columns of the table give the parameter d and 
the respective classical upper bound. The next two columns 
give the best violation of jj" that we have found using the 
d-dimensional maximally entangled state l^'j") in conjunction 
with, respectively, the subset of symmetric multiport mea- 
surements and arbitrary POVMs. The fifth and sixth 
columns of the table give, respectively, the best lower bound 
(LB) and the best upper bound (UB) on the maximal quan- 
tum violation of Jj" that we were able to find. The highest 
level semidefinite relaxation [^, |56[ ^ that we have used to 
obtain the UB is listed in the last column of the table. 
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<5lhv 


5q;i'.„p(i*^)) s 




LB 


UB 


Level 
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0.6667 


0.7124 


0.7124 


0.7124 


0.7124 


2+ 
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0.4800 


0.5366 


0.5375 


0.5376 


0.5578 


1 


7 


0.3878 


0.4587 


0.4587 


0.4587 


0.4668 


1 


11 


0.3058 


0.3325 


0.3325 


0.3328 






13 


0.2781 


0.2987 


0.2987 


0.2991 







For other values of d with d < 13, we have nonethe- 
less found larger quantum violation of b y combining 
the iterative method described in Ref. |l9[ | in conjunc- 
tion with the lower bound (LB) algorithm introduced in 
Ref. Specifically, the following steps were repeated 

a number of times to obtain a non-trivial lower bound on 
the maximal quantum violation of 

(1) Generate alternatively between (i) |\I'J) and (ii) a 
random bipartite pure entangled state in (g)C^; 

(2) Find the best violation and hence the optimal mea- 
surements (for the generated state) using the LB 
algorithm |8|; 

(3) Construct the Bell operator B from the measure- 
ment operators obtained in (2) and determine the 
best violation possible (for these measurements) by 
computing the largest eigenvalue of B: 

(4) Find the best violation and hence the optimal mea- 
surements for the eigenstate [ p9[ corresponding to 
the largest eigenvalue obtained in (3); 
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(5) Repeat steps (3) and (4) until the best violation 
found converges to the desired numerical precision. 

Explicit value of these lower bounds for d < 13 can be 
found in Table III (see Appendix |c| for the quantum 
strategies that realize some of these violations). From 
the table, it is clear that /j" and hence [cf. Eq (H)] 
for d < 13 can be violated by quantum mechanics us- 
ing only projective measurements and l^'J)- This is to 
be contrasted with the results presented by Ji et al. |Q 
where they did not find any legitimate quantum violation 
of their inequalities for = 7, 11 and 13 using mutually 
unbiased measurements. For smaller values of d, it is 
worth noting that our best quantum violation for d = 3 
agrees with that presented in Ref. ||3^; but for d = 5, 
the best quantum violation that we found is about 2.6% 
stronger than the one presented thereof. 

Also included in Table III are upper bounds on the 
maximal violation of I'^ obtained using the semidefinite 
relaxation techniques discussed in Refs. js^, and 
Of particular significance is the upper bound presented 
for d = 3, obtained by considering all the second level 
operators in the hierarchy introduced in Refs. fS^, 15?" 

plus all operators of the form AiB°^B° 

A2B°^By. This upper bound matches exactly the best 
lower bound known, thereby proving that the maximal 
quantum violation of can be obtained using mutually 
unbiased measurements [ p8[ . However, we do not know 
whether the upper bound obtained for d = 5 and d = 7 
can be saturated using quantum strategies. 



CONCLUSION 



is facet-inducing if and only if is facet-inducing. 
We have also investigated the quantum violations 
of I'^ and compared them against those established in 
Ref. [^8| . In particular, we prove that once we are 
equipped with , the corresponding correlation analogue 
proposed by Ji et al. is no longer relevant (in terms 
of detecting an entangled states). On the other hand, we 
do not know if jj" are still relevant once we are equipped 
with the class of CGLMP inequalities . 

In contrast with most other known tight, multiple- 
outcome Bell inequalities |23, |23, f 



11 



and J+ 



are apparently not always violated by a non-maximally 
entangled state. In particular, among the facet-inducing 
inequalities investigated, the best quantum violation that 
we have found for d = 7 is actually due to a maximally 
entangled state. 
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APPENDIX A: AND IJ 



FOR NON-PRIME d 



For completeness, we will also include our results of 
computational investigation in relation to the maximum 
classical value and the tightness of and I'^ for non- 
prime value of d with d < 12 in the following table. 



In this paper, we have reexamined the class of bi- 
partite, d-setting, d-outcome Bell correlation inequali- 
ties proposed in Ref. When rewritten in terms of 
joint probabilities, we show that the corresponding Bell 
inequalities for probabilities naturally generalize the clas- 
sical winning probability of the CHSH game introduced 
in Ref. Q. These Bell inequalities for probabilities, de- 
noted by , are thus also of interest independent of their 
correlation counterpart. 

In establishing these Bell inequalities explicitly, we 
have found that for the more general scenarios of prime 
value d > 5, the authors of Ref. |^ underestimated the 
actual classical upper bounds. Although we could deter- 
mine the actual classical upper bound for d < 13 and 
have provided simple algorithms to estimate them for 
larger values of d, the general problem is left open in the 
present research (the closely related problem for a special 
class of two-outcome Bell inequalities, namely, the XOR 
games is known to be nondeterministic polynomial- 
time hard (NP hard) |§). 

Computationally, we have investigated the tightness 
of the probability inequality and its correlation coun- 
terpart I^^ for d < 13. Our results show that most of 
these inequalities for prime value of d are facet-inducing 
for their respective classical correlation polytopes. How- 
ever, none of these inequalities for non-prime value of d is 
tight (Appendix In this regard, another open problem 
that follows from our observation is whether for each d, 
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38 


12 


1584 


48 


non-tif 


;ht 


17688 


576 


non-tif 


;ht 


60 



APPENDIX B: CLASSICAL STRATEGIES 

Here, we will provide examples of extremal classical 
strategies that realize the values of A^J^xi ^max* and A^^} 
presented in Table |. We will adopt the notation that the 
k-th entry of the vector Oa represents Oa(fc), with the 
exception of Oa(0) which is given as the last entry of the 
vector; likewise for Ob- For d = 7, we have 



Oa 
Ob 



(0 1 2 5 0), 
(0 5 1 6 3 0); 



for d = 11, we have 



o<j = (0 1 10 8 2 5 7 0), 
Of, = (6 7 10 1 5 9 2 3 0); 
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for d = 13, we have 



Table II] can be achieved using l^*^) together with: 



Oa = (0 1 6 8 11 6 4 5 9 0), 


^1 = 


1 

14 


(8, 


7, 


10, 


2, 


10, 


12) 


Ob = (12 6 7 6 7 11 11 10 4 2 0). 




1 

14 


(12, 


3, 


0, 


2, 


8, 


10) 






1 

14 


(0, 


9, 


12, 


8, 


10, 


10) 






1 

14 


(0, 


11, 


4, 


6, 


2, 


12) 


APPENDIX C: QUANTUM STRATEGIES 




1 

14 


(12, 


9, 


4, 


10, 


12, 


2), 






1 

14 


(8, 


3, 


12, 


6, 


12, 


8), 


In this Appendix, we will provide the Schmidt coef- 


^0 = 


1 

14 


(2, 


7, 


0, 


8, 


2, 


2). 



ficients c^''^ 



of the quantum 



state that gives rise to the best violation that we have 
found in Sec. The corresponding quantum states can 
then be written explicitly through their Schmidt decom- 
position, - v-'^-i.Mi 

the phase factors needed to achieve 5, 



J2i=o '^i'\''-A)\'i'B)- We will also provide 

.(1*1)), the 



QM, 

best violation of that we were able to find using the 
d-dimensional maximally entangled state in conjunction 
with the measurements facilitated by a symmetric multi- 
port beam splitter. For this kind of measurements, Alice's 
and Bob's POVM element can be written, respectively, as 
= (;7jc)t n„„ C/^" and B"^ = (U^^)^ n,, U^^ where 
IIo„ = |oa)(oa| and the unitary operators are given by 

= E 4^ ° ''^('^+'^'°)|A:^)(^^|, (Cla) 

u^'=J2^c^'<'^+'^:')\k^){U (cib) 

Note that for each Sq and Sb, we can — without loss of 
generality — always perform the transformations </Jj''° — > 
V?^" —Vo" 1 4>i'' ~^ 4'i^~'t'o' to make (^q" = 0q'' = for all Sa 
and Sb while leaving all the joint probabilities p°_^^'' {sa,Sb) 
unchanged. This is the convention that we are going to 
adopt. In practice, the best multiport measurements that 
we have found are those such that the phases for Alice's 
and Bob's measurements are equal. In what follows, we 
will thus only provide the non-trivial phase factors 



-1} ' 



(C2) 



with the understanding that (p"'^ = ip'^'^ for all Sa- 
Explicitly, for d = 5, the optimal state is 



.(5) 



(0.45367, 0.45049, 0.44898, 0.44378, 0.43899). 



The corresponding measurements that give rise to the 
best violation of arc non-degenerate, and consist only 
of rank-1 projectors. On the other hand, the phase factors 

needed to achieve 5qm. „p (1^*5^)) arc found to be 



(0.92207, 0.65271, 

(0.34126, 0.88988, 

(0.94381, 0.47445, 

(0.14166, 0.23868, 

(0.96047, 0.45749, 



0.79652, 0.22729), 

0.23114, 0.36557), 

0.19652, 0.26924), 

0.79915, 0.05082), 

0.99915, 0.84833). 



For d — 11, and 13, it is expedient to decompose 
the optimal phase factors as 1^9*° = Lp"^ + ip^ , where (/j^ 
and are themselves vectors with {d — 1) entries [of. 
Eq. (p3)]. Explicitly, we have, for d = 11, 





T2 (0, 


0, 


10, 


16, 


16, 


16, 


5, 


0, 


21, 


8), 


•Pd = 


^(20, 


18, 


2, 


2, 


20, 


10, 


5, 


8, 


13, 


8), 


Vd = 


^(16, 


10, 


10, 


2, 


14, 


14, 


13, 


0, 


9, 


10), 




^(10, 


20, 


12, 


16, 


20, 


6, 


7, 


20, 


9, 


14), 






4, 


8, 


0, 


16, 


8, 


9, 


2, 


13, 


20), 


v% = 


22 y ' 


6, 


20, 


20, 


2, 


20, 


19, 


12, 


21, 


6), 


fD = 


22 


4, 


4, 


10, 


0, 


20, 


15, 


6, 


11, 


16), 


Vd - 


^(10, 


20, 


4, 


14, 


10, 


8, 


19, 


6, 


5, 


6), 


^% 


^(16, 


10, 


20, 


10, 


10, 


6, 


9, 


12, 


3, 


20), 




^(20, 


18, 


8, 


20, 


0, 


14, 


7, 


2, 


5, 


14), 




12 (0: 


0, 


12, 


0, 


2, 


10, 


13, 


20, 


11, 


10), 



p, = ^(0.74797,0.65473,0.62522,0.73621,0.82604, 
0.02359, 0.36323, 0.92062, 0.82621, 0.36885), 
and for d = 13, 
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^ (1 

26 ^ ' 

26 (^' 

26 ('^' 

^ (1 

26 ' 

26 (23, 

— fl7 

26 \ ' 

26 (^' 

^(25, 

26 (1"^' 
26 (25, 
26 (^' 

^(17, 
^(23, 



21, 1, 

11, 9, 

23, 11, 

5, 7, 

9, 23, 

9, 7, 

5, 11, 

23, 9, 

11, 1, 

21, 13, 

1, 19, 

3, 19, 

1, 13, 



25, 3, 5, 

15, 19, 25, 

23, 25, 7, 

23, 21, 3, 

15, 7, 13, 

25, 9, 11, 

I, 1, 23, 
21, 9, 23, 
7, 7, 11, 

II, 21, 13, 
7, 25, 3, 

21, 19, 7, 

1, 3, 25, 



2, 18, 
4, 18, 
18, 2, 
18, 22, 
4, 0, 
2, 14, 
12, 12, 
8, 20, 
16, 12, 
10, 14, 
16, 0, 
8, 22, 
12, 2, 



7, 3, 

23, 13, 

21, 3, 

I, 25, 
15, 1, 

II, 9, 
15, 23, 

1, 17, 

21, 17, 

23, 23, 

7, 9, 

25, 1, 

25, 25, 



9, 10), 

7, 10), 

9, 12), 

15, 16), 

25, 22), 

13, 4), 

5, 14), 

1, 0), 

1, 14), 

5, 4), 

13, 22), 

25, 16), 

15, 12), 



In the case of d = 7, the best violation given in 



tp^ = —(0.26436,0.24549,0.26436,0.93681,0.24549, 
26 

0.24549, 0.84021, 0.84020, 0.26436, 0.93681, 
0.84020,0.93681). 

In these two cases, the quantum states that give rise 
to the best violation that we have found can be obtained 
by determining the eigenvector corresponding to the max- 
imal eigenvalue of the respective Bell operator B. Explic- 
itly, these quantum states admit the following Schmidt 
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coefficients: 

c^"' = (0.31463, 0.31456, 0.31352, 0.30525, 0.30462, 

0.30432, 0.30116, 0.29086, 0.29048, 0.28915, 
0.28618) 

c^"' = (0.29189, 0.29189, 0.29189, 0.27790, 0.27790, 

0.27790, 0.27502, 0.27502, 0.27502, 0.27329, 
0.27329, 0.27329, 0.24849). 
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